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By H. Mutt ray 



SUMMARY 

On the "basis of a systematic investigation of G-ottin- 
gen wind-tunnel data on Joukowski type and related air- 
foils, it is shown in what manner the profile drag coeffi- 
cient is dependent on the lift coefficient. It is found 

that, up to camher parameters j- = 0.25 and thickness 

rl 

parameters = 0.4, the profile drag coefficient is rep- 
resented as a cubical parabola whose apex for cambered 
airfoils with c a values lies above zero. With given 

camber these peak c a values are in approximately linear 
relationship with the thickness parameters. The individu- 
al factors for the construction of the profile drag polars 
are given. They afford a more accurate calculation of the 
performance coefficients of airplane designs than other- 
wise attainable with the conventional assumption of con- 
stant profile drag coefficient. 



1. INTRODUCTION 



The assumption of constant parasite drag coefficient 
in the performance calculation of airplane designs (refer- 
ence 1) is particularly inadequate when the c a values 
are high, because a portion of this drag, the so-called 
"profile drag",** upon approaching c a = 1.0, rises consid- 
erably with increasing c a values. As a result the aero- 



*Uber die Abhangigkeit . des Prof ilwider standes vom Auftrieb 
bei Joukowsky- und joukowsky-ahnlichen Profilen." Luft- 
f ahrtf or schung, December 5, 1934, pp. 165-173. 
**The presumption of constant coefficient of residual drag 
(body drag and protruding parts) is preserved. 
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dynamic characteristics, such as criterion of clim"b and 
lift/drag ratio, approach for flight with "best climb and 
ceiling factor rather the values for flight with >egt fine- 
fineness ratio, and consequently do not reach the fig- 
ures arrived, at . with an assumedly constant parasite drag 
coefficient. 

However, since this assumption results in a compara- 
tively convenient method of calculation, it was attempted 
to extend these methods at least to a group of such wing 
sections for which the assumption of constant parasite 
drag coefficient does no longer hold. To this end, we in- 
troduced in Schrenk' s report (reference 2) a second law 
for the quantity of the profile drag coefficient as func* 
tion of c a , to which the representation of the profile 

drag polars as a parabola symmetrical to the c w axis cor- 
responded. Then the profile drag is: 

c 

w prof ile 



In this manner the prof ile • drag coef fie ient con- 
sists of a part independent of c a and a part dependent 

on c a squared. She latter was then considered part. of 

the induced, drag, because' its dependence on c a is as to 

the square also. As a result,' it was possible tointro-- 

duce a substitute span with which the previously derived 

formulas could be analyzed. 

2. THE PROFILE DRAG- LAW 



Now it is readily seen that the law for the profile 
drag contained in formula (l) cannot be carried into ef- 
fect except in very few cases. To begin with, the occur- 
rence- of minimum profile drag at c a = 0 is obviously 
precluded as far as cambered airfoils are concerned. One 
need only think. of the conditions existing on greatly cam- 
bered prof iles to realize this. With such profiles the 
flow, as is known, has already broken down on the bottom 
camber when c a = 0 , so that the polar reveals a break •. 

above c a = 0. But owing to the rounded-off nose* it must 
be preceded by a gradual breakdown of the flow when - com- 

*The validity, is ■ limited to profiles with rounded-off nose. 



' w prof ile min 



+ constant c n 2 



(1) 
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ing from high c a values - approaching the c a value of 
the lower "break (fig. 2) a Thus the conditions, are simi- 
lar to those encountered upon approaching the upper hurt** 
ling point. • Consequently, while" there may exist a para- 
bolic law for the. dependence of profile drag on* lift coef- 
ficient, the apex of the parabola. for symmetrical wing 
sections, :however, can only lie with ca =..'0. Dropping, 
in addition, the assumption, of .squared parabola, the. gen- 
eral term f.or the. prof ile drag then, is: 

%rof" C vr p rof ffiin +C (c a- C a sy /f% rofmi /^ : <?.> 

Hereby, c & = c a value at which the apex of the 

s y m parabola lies 

n = exponent of parabola 
Ac a = c a ~ c a sym 
c = constant 

c w = minimum profile drag coefficient ly~ 

P rof min ing at c~ 

& a sym 

The formula contains only the amount of the power. 

It is now necessary to check whether the given term 
actually agrees with the given conditions, at least of 
normal wing sections, and for the most important flight 
range. 

3. INVESTIGATION OF: G&TTINGEN AIRFOIL DATA 
REGARDING THE. VALIDITY OF THE CITED PROFILE DRAG LAW 



The writer, in collaboration with L. Maxen and F. 
Freytag, undertook the investigation of the validity of 
the law contained in formula (2) on a series of wing sec- 
tions. In order to he able to embody the chosen airfoils 
in a system, the analysis was limited to Joukowski type 
and related airfoils, shown in figure 1 and tahle I. 
This work thus forms a complement to Schrenk* s report (re 
erence 3) and is equally confined to wind-tunnel tests. 
Even the method of defining the profile drag polars by 
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forming differences from the. total drag obtained' by weigh- 
ing, and the induced drag known- frp fa; theory was. preserved 
in .principle . One dif ference, = however. , which is of prima- 
ry .importance for the. size of- the. prof lie., drag coefficient 
and for the establi shment: of . an experimental profile drag., 
law, is; that the calculation of t.ho induced, drag proceeds 
from ..a- rectangular rather than- an e-llipt ical lift distri-. 

•a . 

bu't'ion. JPor equal aspect ratio ~ =1:5,* • there is, as 

known (reference 4) , a 1.04 times induced drag with rec- 
tangular,' as compared to elliptical, .lift distribution. 
Thus at higher c a • the values for the profile' drag would 

"be lower than otherwise generally assumed. For this rea-s 
son we added an increment of 2 percent to the induced 
drag, and the calculation was actually made with 1.06 
times the induced drag obtainable with elliptical lift 
distribution. The reason for this increment was the fact 
that, according to pressure-distribution measurements on 
rectangular wings, the edges reveal high lift peaks, i.e., 
deviations from the theoretical rectangular lift distribu- 
tion t 

; The thus obtained profile drag polars were then di- 
vided in a quota c w not dependent on ca and a 

P r01 min 

quota.. ^ c Wp ro f dependent on c a (fig* 2). The next 

step was the determination of c a . Hereby it was ob~ 
■ a sym 

served that the quota dependent on c a at high Ac a val- 
ues were agreeable as closely as possible because the in- 
strumental errors during the division of the profile drag 
are comparatively smallest in these parts of the polars. 
The third step consisted of plotting on logarithmic scale 
the quota Ac w Qf against Ac a . 



*The fundamental measurements for the tests were made on 
normal wings of -~ = 1:5 aspect ratio. 
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TABLE I. Investigated Airfoils Similar to 
Joukowski Type Airfoils 



Airfoil 
No . 


•p 

■L 

•7 • ■' ' 


d 

0 


. 410 


.0.0 


0.145 


494 


.1 


.0355 


495 


■ *• • .1 


.0 55 


400 


.105 • 


.0583 


496 


.1 


.084 


. 497 


.1 


. 113 


430 


.103 


.113 


498 


.1 


.142 


413 


.105 


.145 


500 


.12 


.055 


502 


.12 


.142 


656 


.2 


.148 


657 ' 


.2 


.148 


523 


.198 


.165 


652 


.1965 ' 


.1665 



4. RESULT OF INVESTIGATION 
a) General result 



The general result of the investigation of Joukowski 
type and related airfoils, which may "be denoted as the 
"profile drag law", may he expressed as follows: 

The profile drag polars for. airfoils with thickness 

parameters ~ ^ 0.4 and camber parameters — = 0.25 may 
t 1 

be replaced with cubical parabolas. The apex of the pa- 

" ered airfoils 1 

b) Magnitude of c 
The obtained c 



rabola for cambered airfoils lies above c a = 0. 



w prof min 

7„.,~. . values for Joukowski and re- 
"proi min 

lated airfoils are shown in figures 3 and 4 plotted against 
i and £. During the determination of the values it was 

» 

found that the profile drag polars occasionally yielded 
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only the values- up- ' to. ■ 0 ..2.5 ..: l3Bicav.se: "by. .definition 

c w . is to "oe the minimum value of the profile drag, 

P rof min , r s ...... -• ' 

which lies at c a .. . .For- .great ly •-cambered airfoils " 

>0.25^ it was found, however, that the lower "break of 

the polar lies above c a , as ' ; a result of which'it was 

impossible to det ermine the magnitude of c_ . 

■■ ♦ - ... , a sym 

An examination of the curved manifests a linear rela- 
tion between c_ and the "airfoil thickness, valid 
.. ; w profmin ; 

up to about ".c™ - = 0.018, respectively, comprising 

P r °f min 

the airfoils'- up to ~ - 0,2 and > ~ = 0.25-. Beyond ^ = 

0.25 the c w _ values rise considerably. The slope 

w .prof m i n " 

• - ' " ' / d \ 

of the straight line c w ; '•'=; f f - ] is unaffected by 

proi m i n .. \l ; 

~. The straight lines may be -expressed with: 

v 

c~ „ - . = f = cvr " , /d \ + 0.044 ~ (3) 

w profmin I w prof m ± n f J = 0 .0 ) I 

The value of c w „ /a ~\ is not 0.0046 for non- 



;w prof min ^-0.oy .. 



cambered airfoils, as the extrapolation of curve 

r _ ff . ., v J\ f (^) reveals. It would be 0.00 55, ac- 
prof m in(j- f O.0 l - VI J 

cording to measurements, on a flat plate at R = 4 X 10 5 , 

the same at, which the airfoil measurements were made. Of 
course-, this figure is valid only under the premise of 
turbulence; of boundary layer over the whole plate depth, 
an assumption which is perhaps not altogether complied 
with'- in the airfoil measurement's. 

The plotting of the c v „ values against -p 

"prof min ~ «- 

(fig. 4) discloses a parabolic relationship of the form 



c w „ /f n \ + const ft^\ (4) 
w prof min |5-= o.O > KiJ 



' w prof min " w prof m i n ^-=o.O^ ' ™"" {J J 
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A logarithmic diagram of the. values gave a parahola expo- 
nent . m <A 3, a figure valid. up to ahout- -^ =-0.225. Be- 
yond it, m increases consistently. 

The values of the constants are found from table II. 



TABLE II. Constants of Formulas 
(4) and. (5) versus ~ 



d 


constant 


0.0-0.2 


. ~ 0.45 


' .25 


~. 50 


.30 


~. 60 


.40 


' ~.90 



as 



Now Co . versus r* and -r may be approximatec 

P roi min . * 



c w 0.0046 + 0.044 r + const... 



3 

(0 



(5) 



The formula is valid, as stated before, up to j = 



0.25 and ^ = 0.225. Putting- the constant for this range 
at 0,48, we obtain 



c w = 0.0046 + 0.044 + 0.48 

w prof mill I 

as approximation formula. 



(!) 



(5«) 



c) The Magnitude of. Ac w . 



prof 



Ac w f is, as previously explained, a function of 



Aca = c a - c. 



sym* 



for which reason the determination 



of the magnitude of c a sym formed the second step in the 

definition of the profile drag. However, in order to com- 
plete the discussion of the quotas of the profile drag 
first, it is deemed best to revert to the magnitude of 



^° w prof an<i hereby assume the magnitude of c a 
known. 



sym 



to be 
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Figures 5a to 5f and 6a' to" 6f show as' examples the 
logarithmic ' diagrams of the Ac w " • values ver.sus Ac a for .. 
two series of regular Joukowski airfoils:, once for con- .. 

stant thickness (%■ = 0.1) and variable camher (% = 0.0 

- • - v ^ . / ... ..... ... /f A 1 

to 0.25 j and then for constant camher = 0.1 J and 

variable thickness (r ~ 0.0 to 0.4'-^. ■ -A straight line 

with 1:3 slope is drawn through the points. Naturally 
this method of plotting gives the points a more or less 
pronounced scatter because the: instrumental errors now en- 
ter only in the Ac a . The points for high Ac a values 
have been particularly observed when drawing the straight 
lines. They lie quite well on the straight lines. Strict- 
ly viewed, the slope of the straight. lines naturally fluc- 
tuates also, but the coordinated constant (intersection of 
straight line with the parallel to axis Ac w at Ac a = 
1,0) of the chosen slope, affords a certain equalizat ion* 

The remaining figures may equally well be obtained 
with straight lines of 1:3 slope. All measurements, of- 

course, are subject to the limitation that |- must not 
exceed 0.25. 

The c constants of formula (2) defined as straight 
lines in logarithmic • diagram, are shown in figures- 7 and. 
8 for pure and modified Joukowski wing sections versus £ 

j "-<<•■■' 
and There is a pronounced scatter which makes it dif- 

ficult to obtain curves with' Very many test points, al- 
though it was quite satisfactory in isolated cases as, for 

example, in the plotting of c versus j- for ^ = 0.1, 

0.15, and 0.3. The other curves were drawn accordingly, 
without regard to any marked, or perhaps, even one-sided 
scattering of the corresponding test points. 

The plotting of the c constants against 4. proxi- 

mated hyperbolas, and against |- parabolas. Because of 

the pronounced scatter, no formal evaluation of the curves 
was attempted. 
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d. Magnitude of c a sym 

The determination of this value is of particular sig- 
nificance insofar as it is possible theoretically to de- 
termine the location for infinitely thin airfoils, and it 
was found that this location was, as anticipated, in ac- 
cord with that obtained "by experiment. Further, it was 

established that, as regards c a . there exists a sim- 

»sym 

pie corelation for the airfoils of finite thickness. 

Let us first examine the experimental values (figs. 

9 and 10). Figure 9 shows that the c a values for 

°sym 

airfoils of equal camber plotted against |s may be ap- 
proximately connected by straight lines. These lines are 
so much steeper as the camber is greater. The origin of 
the straight line in the ordinate axis is coincident with 
the point of c a (obtained by extrapolation) for 4 = 

0.0. sym * 1 

The extrapolation was made as follows: From the 

c a ., and c fl diagrams -plotted against a for 

**theor a exper 

the Joukowski airfoils, we took the conversion factors 

independent of a, and plotted them against the 
c atheor 

camber for j = 0.05 to 0.25 (fig. 11). The result was 

a set of parallelly displaced curves. It was assumed that 
this parallel shifting applied very likely to the curve 
cL 

for — = 0 also. For obvious reasons the origin of the 

curve for f = 0 had to be in point --$££.2*- = i. with 
I ^ a theor 

these values of the quotients the theoretically obtainable 
values of c a sym for f" = 0 then gave the corresponding 
experimental figure. 

No definite statement may be made about the curves 
in figure 9 within range of |. > 0.4, because of the 

lack of experimental data on airfoils of greater thick- 
ness. Consequently, the extrapolation of the straight 

line over f = 0.4 in figure 9, is not absolutely cer- 
tain. 1 



a exper 
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' Figure. 10 "gives' . c a ' ; ye* sua According, to it. 

the" rise at ' fir'st'- is : subsequently' followed' /by" a drop in the 
values within range of the! tests for. . ^. >,'6,25. . Hence the 

r e-currence ' ©f' • Co ■ •- = 0 for the- same profile' thickness,. ■ 

a sym ' " 

once for. ™ =• 0 and again for: • |- > 0 . ' •. 

' IText.. it -was attempted to . esta.bli.sii one .particular . . 

point on the airfoil to which the pertinent value* of 

c n . could lie. somewhat- coordinated.: ■■. - 
a sym • - 

An. 'excellent point for'" cambered 'airfoil's' of zero 
thickness ia the leading; edge . .' In order . to obtain, the 
minimum /prof iie . drag the "forward, stagnation point .must " ey- . 
identlylie in' the leading edge' because , .- select ihg a point, 
in the vicinity of the leading edge, the circulation of' " 
the leading edge, respectively the ensuing great veloci- 
ties, would undoubtedly" involve ; a much higher profile drag 
than' with the so-called " shock-free ' entry" of the flow,.. 
Adding . further the ..trailing .edge in the- known manner as - : 
rear stagnation point,' the' angle of attack of the airfoil 
and . the corresponding ,. c a and Ca s ym value (reference's) 

is known according to the equally assumed as known re- 
sults of the' conformai: transformation. 

Passing to profiles of finite thickness, one might be 
tempted to' assume that the forward stagnation point to 
which 6 a sy in - is coordinat ed , ■ is coincident with the inter- 
section of. the profile contour and the extension of the 
skeletbn line of. the thick airfoil (profile with zero 
thickness - and.- equal camber)*. But the design figure (fig. 
12) shows that' a substantially higher theoretical angle 
of attack belongs to this point (point' 2); that is, a much 
higher lift also. On the other hand, according to figure 

9, the e a „„„ value decreases as the profile thickness 
' a sym c 

increases. This signifies" that the forward stagnation 
point of the cited intersection of skeleton line and pro- 
file, .contour (point 2 in- fig. 12) moves forward (in direc- 
tion of the arrow)., Such .a. po int . of • the profile is the 
intersection of profile contour and axis of. abscissa 
(point 3 -in fig. 12) i We coordinated the same theoretical 
angle of attack to this point independent of y, although 

the corresponding lift coefficient still continues to rise, 
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albeit slightly, as a result of the increasing thickness 
of the profile. The pertinent c a values for. this point 
are included in figure 9. 

As previously stated, there is no agreement between 
experimental and theoretical values for cttheor = const. 

For this reason the corresponding curves were not plotted 
in figure 10, although it contains the theoretical curve 
for ^ = 0.0. 

Graphical determination of the actual position of the 
forward stagnation points by means of the cc sym values 

coordinated to the experimentally defined c a values; 

a syui 

that is, on the "generating circles" of the profiles as 
well as on the profile contours themselves, discloses 
these points in the design figures to again lie approxi- 
mately on straight lines (figs. 13 and 14), The test 
points in the Z plane, that is, in the representation of 
the "generating circles" are best replaced with a straight 
line, whereas this is hardly possible for greater |- in 

the representation of the Joukowski airfoils. One obtains 
instead slightly cambered curves in the Z plane. 

The straight lines containing the forward stagnation 
points in the representation of "generating circles" are 
steeper as the profile camber is greater. It ties up with 
the relationship existing between figure 9 and figures 13 
and 14. But the advantage of the latter over figure 9 is 
their lucidity as well as the fact that the scattered test 
points may be closely replaced by straight lines through 
this method of representation. Only by computing the 
° a sym val ' ae s corresponding to these straight lines were 

the slightly cambered experimental curves of figure 9 and 
the experimental curves of figure 10 established. 

How the theoretical values may be obtained from the 
experimental for any thickness parameter with the quo- 

. ° a exper 

tients — of figure 11, in the same manner as the 

a theor 

Ca sym eX p for T ~ °«°» 111(3 curves obtained in this fash- 
ion are shown as dashed lines in figures 9, 10, 13, and 
14 for the profiles with maximum £ and ^. 
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Figure .15: finally shows, the pitch of the straight 
lines connecting the forward; stagnation pointts of figures 
13' and 14' plotted against the camper. The values diverge 
comparatively little from a parabola. 



5. APPLICATION OF THE RESULTS TO THE PERFORMANCE 
• CALCULATION OF AIRPLANE DESIGNS 



In the ' numerical : example hereafter, we show what dif- 
ferences occur in the magnitude • of the most essential 
aerodynamic performance factors when the caloulat ion is 
made with cori'stant profile drag coefficient and when ef- 
fected' with-prof ile drag coefficient dependent or. lift 
coef f icient . ' : As example, we use a glider with non-twisted 
elliptical wings "of 1:15 aspect ratio and constant pro- 
file; wing warping and profile change of outer wings shall 
be disregarded for reasons of simplicity. The chosen pro- 
f ile. is a Joukowski airfoil with a thickness paramoter 

— ='0»i'25, camber parameter y = 0.15, and wing span = 15 
t » 

motors. 



c w „ 0.012 
w prof min . 

c = .0325 



C a _ — = 



sym 



Figure 3 yields 

Figure . 11 " 
Figure 13 » . 
As residual drag coefficient, we take 
c Wr = 0 .0025 

'Smin 



so that c Wc . . = .0145 



therefore, f w = .0145 X 15 = 0.2175 m 3 

s min 

Whence with point-by-point evaluation, the lift coeffi- 
! a : 



T/ith best criterion of climb at c a + = 1.05 (i.4325) 
it 



' a st 

" lift/drag ratio at c a£ ^ .80 (0.826) 
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The criterion of climb — ~r- = s 

c a 3 3 



With best criterion of clime at ;s gt = 3 7" ' 375 3 9 ' 375 



" lift/drag ratio at 



1 .1 



25.475 25.84 



c. 



The lift/drag ratio -- = S 



With best criterion of climb at ' £ st = 26.8 (24.55) 

" 11 lift/drag ratio at e £ = 28.5 (28.45) 
The conversion factors are:' 

c a D + 9 st 

— -- = 1.313 (1.734); - — = 0.93 (0.882) 
c a e 4 € 

£st _ 0#941 (0.8625) 

The figures given in parentheses are the quantities 

obtained for constant parasite drag coefficient c w = 

s min 

0.0145 according to Schrenk (reference l) . The "bracketed 
figures for the case of "best lift/drag ratio differ very- 
little from those corresponding to the actual course of the 
polar. But the figures for the best criterion of climb, on 
the other hand, reveal a perceptible difference, especially 
the c a values. The c a value of the best climb and ceil- 
ing factor ( c a s t = 1*4325) computed with constant profile 

drag coefficient, lies about 36.5 percent too high; that 
is, so high as to actually be quite close to the value of 

c a = 1.5 of the airfoil. 

"max 



Translation by J. Vanier, 
national Advisory Committee 
for Aeronautics. 
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Figure 3.- Minimum profile drag versus thickness parameter (d/i) 
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Figure 5 c. 
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Figure 5 b. 
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Figure 5 d. 
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Figure 5e. 
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Figure 5f« 



Figure 5(a-f)«- Logarithmic plotting of profile drag quotas dependent on 
lift of Jouknweky airfoil* with d/Z = 0.1 
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Figure 6 a. 
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Figure 6 h. 
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Figure 6 c. 
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Figure 6d. 
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Figure 6e. 

Figure 6(a-f).- Logarithmic plotting of profile drag quotas dependent on 
lift of JotfcowBky airfoil* with f/l = 0.1 
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7.- Relationship of c constant of 
formula (2) to d/Z . 
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Figure 8.- Relationship of c 
constant of 
formula (2) to f/z . 
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Mgure 9.- Relation of lift coefficient Figure 10.- Relation of lift coefficient 
for minimum profile drag for minimum profile drag 
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Figure 12.- Representation of travel of 

forward stagnation points 
corresponding to minimum profile drag 
with increasing profile thickness. 




Figure 13.- Position of forward 

stagnation points for 
minimum profile drag on the 
"generating circles" and on the 
Joukowsky airfoils themselves 
(f/Z » 0.1). 
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Figure 11.- Conversion factor Cagxp/^theoy 
of Joukowsky airfoils versus f/l . 



Figure 15.- 

Angle 6 
of figs. 13 
and 14 
versus f/l . 
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Figure 14.- Position of forward 

stagnation points for 
minimum profile drag on the 
"generating circles" and on the 
Joukowsky airfoils themselves 
(f/l » 0.2). 



